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Abstract 

We consider a reaction-diffusion equation in a cellular flow. We prove that in the strong 
flow regime there are two possible scenario for the initial data that is compactly supported 
and the size of the support is large enough. If the flow cells are large compared to the 
reaction length scale, propagating fronts will always form. For the small cell size, any finitely 
supported initial data will be quenched by a sufficiently strong flow. We estimate that the 
flow amplitude required to quench the initial data of support L is A > CLq\ii(Lo)- The 
essence of the problem is the question about the decay of the L°° norm of a solution to the 
advection-diffusion equation, and the relation between this rate of decay and the properties 
of the Hamiltonian system generated by the two-dimensional incompressible fluid flow. 

1 Introduction 

It has been well understood since the classical work of G.I. Taylor that the presence of a fluid flow 
may greatly increase the mixing properties of diffusion. This phenomenon is known as "eddy 
diffusivity" or "enhanced diffusion". The mathematical approach to the problem is usually via 
the homogenization techniques that concentrate on the long time-large scale behavior: see 
for a recent extensive review. This approach is appropriate when there are no other time scales in 
the problem so that one may wait as long as needed for the mixing effects to become prominent. 

Recently there has been a lot of interest in the effect of flows on the qualitative and quantita- 
tive behavior of solutions of reaction-diffusion equations. Intuitively, there may be two opposite 
effects of the additional mixing by the flow: on one hand, it may increase the spreading rate of 
the chemical reaction (the "wind spreading the fire" effect), or it may extinguish the reaction 
(the "try to light the campfire in a wind" effect). The first effect is related to the behavior of 
front-like solutions, and has been extensively studied recently: traveling fronts have been shown 
to exist in various flows El 13 EI3 EH ESI EH , and flows have been shown to speed-up the front 
propagation due to the improved mixing EJ El 1213 EZ| j see [U EEj for recent reviews of the 
mathematical results in the area. This problem has also attracted a significant attention in the 
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physical literature, we mention El El Ed EEE1 ECU 1221 E3J among the recent papers and refer 
to jHU] as a general reference. The present paper addresses the second phenomenon mentioned 
above: the possibility of flame extinction by a flow. The basic idea is that if the reaction process 
may occur only at the temperatures T above a critical threshold 8q, then mixing by a strong 
flow coupled to diffusion may drop the temperature everywhere below 9q an d hence extinguish 
the flame. However, unlike the usual linear advection-diffusion homogenization problems, one 
may not wait for this to happen beyond the time t c it takes for the chemical reaction to occur - 
the mixing has to happen before this time. The question we address is: "Given a threshold 6q, 
a time t c , and the support Lq of the initial data, can we find a flow amplitude Aq(Lq) so that if 
the flow amplitude A > Aq(Lq) is sufficiently large then sup x T(i = t c , x) < #o?" This problem 
has been first considered in 10 for unidirectional, or shear, flows that have open streamlines. 
Even in this simple situation the answer is non-trivial: in order for quenching to be possible, 
the profile u(y) should not be constant on intervals larger than a prescribed size. The answer 
has been shown to be sharp in 21 . 

In this paper, we study quenching by a different class of flows with a more complex structure: 
incompressible cellular flows. These are flows such that the whole plane R 2 is separated into 
invariant regions bounded by the separatrices of the flow that connect the flow saddle points. 
Many types of instabilities in fluids lead to cellular flows, making them ubiquitous in nature. 
We only mention Rayleigh-Benard instability in heat convection, Taylor vortices in Couette flow 
between rotating cylinders or heat expansion driven Landau-Darrieus instability. The fact that 
the cellular flows have closed streamlines make the effect of advection more subtle. An important 
role in the possibility of quenching is played by a thin boundary layer which forms along the 
separatrices of the flow. Our main results show that the cellular flow is quenching if and only if 
the size of the minimal invariant regions (the flow cells) is smaller than a certain critical size of 
the order of laminar flame length scale. 

The simplest mathematical model that describes a chemical reaction in a fluid is a single 
equation for temperature T of the form 



where the flow u(x) is prescribed. We are interested in the effect of a strong advection, and 
accordingly have written the velocity as a product of an amplitude A and a fixed flow u(x). In 
this paper we consider nonlinearity of the ignition type, that is, we assume that 



The threshold 8q is called the ignition temperature. The last condition in (J2J is just a nor- 
malization. We consider the reaction-diffusion equation JIJ in a two-dimensional strip D = 
{x G R, y G [0, 2irl]} with the periodic boundary conditions at the vertical boundaries: 





T(0,x) = T (x) 



(i) f(T) is Lipschitz continuous on < T < 1, 

(ii) /(l) = 0, 30 o such that f(T) = for T G [0, O ], f(T) > for T G (0 , 1) (2) 

(iii) f{T)<T. 



T(x,y + 2TTl) = T(x,y). 
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The initial data Tq(x) = T(0,x) is assumed to satisfy < Tq(x) < 1. The maximum principle 
implies that then < T(t, x) < 1 for all t > 0. We will say that regions with temperature close 
to one are "hot", and those with temperature close to zero are "cold". 

The problem of extinction and flame propagation in Q with the ignition type nonlinearity 
((21) was first studied by Kanel . 16^ in one dimension and with no advection. Assume for simplicity 
that the initial data are given by a characteristic function: 7b (x) = X[o,L]( x )- Kanel showed that, 
in the absence of fluid motion, there exist two length scales Lq < L\ such that the flame becomes 
extinct for L < Lq, and propagates for L > L\. More precisely, he has shown that there exist 
Lq and L\ such that 

T(t, x) —* as t — > oo uniformly in D if L < Lq (3) 
T(t, x) —* 1 as t — > oo for all (x, y) £ D if L > L\. 

In the absence of advection, the flame extinction is achieved by diffusion alone, given that the 
support of initial data is small compared to the scale of the laminar front width l c = M -1 / 2 . 
However, in many applications quenching is a result of strong wind, intense fluid motion and 
operates on larger scales. Kanel's result was extended to non-zero advection by shear flows by 
Roquejoffre |31| who has shown that (J3J) holds also for u ^ with Lq and L\ depending, in 
particular, on A and u(y) in an uncontrolled way. 

As we have mentioned, the question of the dependence of the strength of advection A which 
is necessary for quenching the initial data of a given size Lq has been recently studied in ^01 and 
|21| in the case of a unidirectional (shear) flow (Au(y),0). Following ^U], we call the flow u(x, y) 
quenching if for every Lq there exists Aq(Lq) such that the solution of Q with the initial data 
of size Lq and advection strength A > Aq(Lq) quenches. It turns out |1() | I21| that the shear flow 
u(y) is quenching if and only if u(y) does not have a plateau of size larger than a certain critical 
threshold (comparable with the length scale M -1 / 2 which characterizes the width of a laminar 
flame). The intuition behind this result is that shear flows are very effective in stretching the 
front and exposing the hot initial data to cool-off effects of diffusion unless there is a long, flat 
part in their profile, where this phenomenon is obviously not present. 

Here we consider for the domain D = R x [— ttI,ttI] with the 27r/-periodic boundary 
conditions in y and decay conditions in x: 

T(t,x,y)=T(t,x,y + 2Ttl), T(t,x,y) — >0 as x — ► ±oo. (4) 

We restrict ourselves to a particular example of a cellular flow u(x, y) that has the form U[(x, y) = 
V hi(x, y), where V^/i = (h y ,—h x ). Here I defines the size of a flow cell, and we take the stream 
function hi to be 

hi(x,y) = I sin j sin y (5) 

whose streamline structure are shown in Fig. 1. 

We will usually omit the index / in notation for m, hi for the fluid flow; it will be clear from 
the context what the scaling is. The initial data To(x,y) is non-negative, and bounded above 
by one: < To(x,y) < 1. 

The first theorem shows that cellular flows with large cells do not have the quenching prop- 
erty. 
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Figure 1: Cellular flow 



Theorem 1.1 Assume that To(x,y) = 1 for (x,y) £ [0, 7rZ] x [0,7tZ]. There exists a critical cell 
size Iq ~ M~ 1 / 2 so that if I > Iq, then under the above assumptions on the advection u, we have 
T(t, x, y) — ► 1 as t — > +oo, uniformly on compact sets, for all 4 

The notation l ~ M" 1 / 2 means that CiM" 1 / 2 < l < CW" 1 / 2 with Ci )2 some positive 
universal constants. The proof of Theorem 11.11 is simple and is based on a construction of an 
explicit time-independent sub-solution. 

Next, we show that if the cell size I is sufficiently small then a sufficiently strong flow will 
quench a flame. More precisely, we have the following result. 

Theorem 1.2 Assume that Tq(x, y) = outside an interval —Lq < x < Lq and < Tq(x, y) < 1 
for all (x,y) G D. There exists a critical cell size Iq ~ M -1 / 2 so that if I < Iq, then there exists 
Aq(Lq) such that we have T(t,x,y) — > as t — > +oo, uniformly in D, for all A > Aq(Lq). For 
large Lq, we have Aq(Lq) < C(/)Lq ln(Lo). 

A formal argument based on the homogenization theory predicts that Aq ~ Lg without the 
factor of In Lq - this follows from the effective diffusivity scaling ~ V~A that was first shown 
formally in [S] and later proved in |X 1 1 1241 128j . The same scaling may be obtained from the 
formal predictions Va ~ A l / A for the front speed Va in a cellular flow P3 121 El E2] ~ this implies 
that the front width is of the order A l / A . Hence one might expect that initial data with the 
support less than the front width Lq < A 1 ^ to be quenched. Therefore, the rigorously proved 
bound of Theorem 11.21 is likely to be sharp up to a logarithmic factor. 
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The particular choice of the stream function © is not important for the proof but it does 
simplify some of the estimates - it is straightforward to generalize our result to other cellular 
flows. However, the proof of Theorem ll.2l does use the periodicity of the flow in an essential way. 
We believe that the quenching property should hold for sufficiently regular non-periodic flows 
with small cells as well. We have preliminary results in this direction using different techniques; 
however, these results give much weaker upper bound for Aq(Lq), and will appear elsewhere. 

At the heart of the proof of Theorem ll.2l is the question about the rate of decay of L°° norm 
(and its dependence on A) of the solution to a passive advection-diffusion equation. Thus, our 
main object of study is a natural question about the effects of a combination of two fundamental 
and separately well-understood processes: advection by a fixed incompressible flow and diffusion. 
Yet their interaction is well known to produce subtle phenomena. The issues we study are directly 
related to the work of Freidlin and Wentzell |13] 1141 115] on the random perturbations of the 
Hamiltonian systems. They show that in the limit of large A the process converges to a diffusion 
on the Reeb graph of the background Hamiltonian. The relation to that work is very natural as 
any incompressible flow in two dimensions is a Hamiltonian system (the stream function is the 
Hamiltonian). However, the cellular flow that we consider does not satisfy the assumptions of 
the Freidlin- Wentzell theory, which requires growth of the Hamiltonian at infinity and does not 
allow existence of hetero-clinic orbits. Nevertheless, one may restate the small cell assumption 
in Theorem 11.21 as a requirement that the Reeb graph of the Hamiltonian has a sufficiently 
small diameter in its natural metric. The proof of Theorem 11.21 is based on two observations: 
first, temperature becomes approximately constant on the whole skeleton of separatrices, as the 
skeleton is just one point on the Reeb graph. Because of that, solution inside each cell may 
be split into two parts. One solves the initial value problem with zero data on the boundary. 
Another solves a nearly identical boundary value problem in each cell. The first part decays 
because the cell is small - this is where we use the size restriction. The second one is nearly 
identical on all cells, hence it has to be small in order not to violate the preservation of the 
L 1 -norm. The technical part of the proof is in making this informal scenario rigorous. 

The paper is organized as follows: in Section |2] we prove Theorem 11.11 by constructing an 
appropriate sub-solution and using certain PDE estimates to prove convergence of solution to 
unity. In Section |3] we give the proof of Theorem ll.21 which is more involved. It uses uniform in 
A estimates on the evolution of advection-diffusion equation, a boundary layer argument, and 
probabilistic estimates for an auxiliary cell heating problem. 

Acknowledgment. The research of AF is supported in part by The Centennial Fellowship 
from American Mathematical Society and U.S. National Science Foundation (NSF) grant DMS- 
0306659. AK has been supported in part by NSF grants DMS-0321952 and DMS-0314129. LR 
has been supported in part by NSF grant DMS-0203537 and ONR grant N00014-02-1-0089. 
Both AK and LR acknowledge support by Alfred P. Sloan fellowships. 

2 Absence of quenching by large cells 

We prove in this section Theorem 11.1] that is, we show that cellular flows with sufficiently large 
cells do not have the quenching property. The proof consists of two steps. First, we construct a 
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time-independent sub-solution $(x) to (Q) in a cell C\ = [0, tvI] x [0, ttI]: the function satisfies 

Aui (x) • V$ < A$ + M/($) (6) 

and is 27r/-periodic in y. It is also positive on an open set inside C\ and negative on dC\. We 
normalize $ so that «3> < 1. As To = 1 on C\ by assumption, we have 7o(x) > $(x). Then the 
maximum principle implies that T(t,x) > $(x) for all t > and x G Ci. It follows that T(t, x) 
does not vanish as i — ► +oo. In the second step we show that actually T(t,x) — ► 1. We begin 
with the construction of the sub-solution ^(x). First, we rescale equation Q by x — > /x so that 
a sub-solution in the rescaled coordinates should satisfy 

4«(x) • V$ < iA$ + M/($). (7) 



Lemma 2.1 // / is sufficiently large then there exists a C 1 function <&(x) that is constant on 
the streamlines of the flow it(x) and satisfies {71) /or x G Ci = [— 7r, 7r] 2 . Moreover, ^(x) < 1 /or 
a// x G Ci, and <&(x) < /or x G <9Ci. 

Proof of Lemma 12.11 We may choose two numbers Q\ and 02 so that 8q < 9\ < 02 < 1 and 

such that the straight line that connects the point (#i,0) to the point (02,f(&2)) lies below the 
graph of /(T). More precisely, that means that the function 

f 0, T < 6*i 

g(T) = \a(T-8 1 ), 0!<T<8 2 (8) 
I /(T), T>0 2 

with a = f{02)/{02 — Qi), satisfies g(T) < f(T). Such modification of f(T) for the construction 
of sub-solutions was first used in |16j . and then in A function $(/i(x)) satisfies if 

where h(x,y) = sin a; sin y is the stream function. Note that the advection term vanishes identi- 
cally for such functions. We use the fact that Ah = — 2h to obtain 

Here Z c = M -1 / 2 is the laminar front width, the length scale associated to the chemical reaction 
strength. Relation © indicates that the ratio l/l c has to be sufficiently large for a sub-solution 
to exist. Note that 

2h(l -h)< |V/i| 2 < 2(1 - h 2 ). 

Indeed, we have 

\v->l,l \|2 -2 2 , 2-2 • 2 , • 2 r> • 2 ■ 2 

| vri{x, y)\ = sm xcos y + cos xsm y = sin i + sm y — 2sm xsm y 
> 2 sin x sin y(l — sin x sin y) = 2h(x,y)(l — h(x,y)) 
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and 



\Vh{x,y)\ 2 < 2 - 2 sin 2 x sin 2 y = 2(1 - h 2 (x,y)). 
Therefore it suffices to construct an increasing function ^(/i) satisfying 



2h d<$> 1 
+ 



dh? 2/i(l -h)dh 2(1 - h?) \k 
which would in turn follow from 



t) #)>o. 



1 d<$> 1 / I 
+ 



Make a change of variables R 



d 2 <f> 

~dh! 1 I - h dh ' 2 V / 
/ 



#) = o. 



kV2 



(1 — h) so that the above becomes 



d 2 $ 1 d$ 



(10) 

z 



ciR 2 RdR 

The center of the cell corresponds now to R = 0, while the boundary /i = becomes i? 
We impose the following "initial data" for (|1U|) : 

If- 

The explicit form (jHJ) of the function g(&) implies that the solution &(R) is given explicitly by 

6 



= ^ + (0 2 - ^)J (fl^/a) , for R < Ri 
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with a as in Q. Here Jo(£) is the Bessel function of order zero, and £1 is its first zero. 
Furthermore, we have 

R-? 

$(i?) = Bln-f , iovRi<R. 
R 



(12) 



The constants B and i?2 are determined by matching the functions (|llj) and ()12|) . and their 
derivatives at R = R\ . Then we get 



B = {e 2 -e l )i 1 \j^ l )\, #2 = 6 



/(#2 



■ cxp 



? 2 -0i)6K(6 



Observe that the function <£(-??) constructed above is negative on the boundary of the cell only 
provided that R2 < ;=, which means that the cell size 



exp 
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has to be sufficiently large for this construction to be applicable. This proves Lemma 12. II □ 

In order to finish the proof of Theorem 11.11 we have to show that T(t,x) — ► 1 as t — > +00 
provided that 7o(x) = 1 on a cell. For such initial data we have 

T (x) > $ (x) = max{$(/j(x)),0} , 

where the function ^(x) is the sub-solution constructed in Lemma 12.11 It follows from the 
parabolic maximum principle that then T(£, x) > $o( x ) f° r an t > 0. Furthermore, we have 
T(t,x) > ty(t,x), where the function ^(£,x) satisfies (^Q) with the initial data <I>o( x )- Note 
that x) > ^olx) f° r a h t > 0. The maximum principle applied to the finite differences 
^fh(t,x) = *$>(t + h, x) — x) implies that *$>(t,x) is a point-wise increasing function of time 
that is bounded above by one. Therefore the point-wise limit ^(x) = limt_ + + 0O ^f(t, x) exists, 
moreover, ^(x) > <I>o( x )! an d ^( x ) satisfies the stationary problem 

Au(x) ■ V\F = A* + Mf(W). (13) 

with the 27r/-periodic boundary conditions in y 

*(t,a;,0) =*(* } x,27rl) 

(here 2N is the number of cells in y direction). We have the following lemma. 

Lemma 2.2 The function *$>(x,y) satisfies the following bound: 

[ |V^(x)| 2 dx+ f /(tt(x))dx < +00, (14) 
Jd Jd 

where D = M x x [—irl,irl]y. 

Proof. The function ^Sf(t, x) satisfies an a priori bound 

r-i 



tJo \J X ^ / d * ~ C ' ' r - CM " 



that may be easily proved as in j^j. Here the constant Co may depend on the flow amplitude A. 
Therefore, there exists a sequence of times t n — * +00 so that 

/(*(tn,x))dx<a,. 



This implies that 

/ /(*(x))dx<C 0) 

JD 

and it remains to obtain the bound on HV^H^ in l|14jl. We multiply (|13j) by ^ and integrate 
in x between —X + £ and X + £ with X large and £ 6 [0, l c ], and in y G R. We get 

4 r [u^X + Cy^iX + Cy^-mi-X + Cy^i-X + Cy^dy 
= f l [W(X + C, + C, V) - *(-X + C, V)*x(-X + C, y)] dy 

J-ttI 

+M / / y(x,y)f(y(x,y))dx- dy \W(x,y)\ 2 dx 

J-nl J-X+C J-nl J-X+C 
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and average this equation in £ £ [0, l c ]. This provides the bound 



1 



c Jo 



nl 

dC I dy 

-nl 



XH _ 2-irl — 

\V^(x,y)\ 2 dx <tiAI\\u\\ 00 + — + M I f(V(x,y))dxdy 



D 



Taking X to infinity, we obtain 



\V^{x,y)\ 2 dx <-kAI\\u\\ 00 + z ^- + M I f(^{x,y))dxdy 



2ttI 



D 



D 



and the bound on V\I/ in (|14j) follows. □ 

Now we are ready to complete the proof of Theorem 11.11 
Proof. Lemma 12.21 implies that there exist two sequences of points x T: 
so that 



-oo and z n 



/nl 
-nl 



W(x n ,y)\ 2 + dy -► as n -► +oo. 



We integrate (|T3|) in y and in x between x n and z n to obtain 



A 



-i 

-irl 
irl 

nl 



[ui(z n ,y)y(z n ,y) - ui{x n ,y)^{x n ,y)] dy 

r— — i r l f Zn — 

\ty x (z n ,y)-ty x (x n ,y)]dy + M dy f(^(x,y))dx. 



nl J x 



We pass to the limit n — > oo in ()16j) . Observe that 



fy x (z n ,y)dy 



-nl 



< VttI 



d \ 1/2 

|2, 



|^(z n ,y)| <iy 



-ttZ 



as n — > +oo 



as follows from (|15j) . and similarly 



7,1 



"S> x (x n ,y)dy 



as n — > +oo. 



Furthermore, since 



Til 



-7T« 



u 1 (x,y)dy = 



for all (|15|) and the Cauchy-Schwartz inequality imply that 



nl 



u 1 (z n ,y)^(z n ,y)dy 



nl 



nl 



nl 



ui(z n ,y) / ^(zn^KUl/^O, asn^+cx). 



Therefore in the limit n — > +oo equation H16|) becomes 

f(V(x,y))dxdy = 



+oo 
(15) 

(16) 



and hence f(^(x,y)) = for all (x, y) G D. However, since ^(x,?/) > &o(x,y), while, on the 
other hand, max ^o(x,y) = 62 > 0o, and the function ^/(x,y) is continuous, we conclude that 
ty(x, y) = 1. This finishes the proof of Theorem ll.il □ 
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3 Quenching by small cells 



In this section we show that the cellular flow with small cells is quenching. The proof proceeds 
in several steps. First, we reduce the problem to a linear advection-diffusion equation. Indeed, 
as fiT) < T we have the following upper bound for T : 

T(t,x) < e M V(i,x). (17) 

The function 4>(t,x) satisfies the advection-diffusion equation 

^ + Au ■ V0 = Acf) (18) 

with the same initial data 0(0, x) = 7b(x) and the 27r/-periodic boundary conditions in y: 
(j>(t,x,y) = cf)(t,x,y + 2irl). Note that if at some time to > we have 4>(to,x) < 6q everywhere, 
then the maximum principle implies that </>(i,x) < 8q an d T satisfies the linear equation (|18|) 
for all t > tQ. Then the conclusion of Theorem 11.21 follows . Hence, the upper bound (|17|) implies 
that it suffices to show 4>(t = M _1 , x, y) < 9' = #o e_1 an d this is what we will do. 

Heuristically, the proof relies on the observation that solution of (|18|) should generally become 
constant along the streamlines of the flow if its amplitude is large. Moreover, the value of the 
solution on the streamlines h = ho very near the boundary in two neighboring cells have to 
be close (as follows from a simple L 2 -bound on X7<fi appearing in Lemma 13.31 below) . However, 
that means that solution should have, roughly speaking, the same profile in each cell. This is 
incompatible with the preservation of the L 1 -norm of <fi unless this function is very small in each 
of the cells which means that solution has to be less than 9q everywhere. The proof follows this 
heuristic outline - the technical difficulty is that we are able to control the uniformity of the 
solution along the streamlines only in a space-time averaged sense. Additional ingredients are 
required to obtain the point-wise control. 



3.1 The Nash inequality lemma 

We will need throughout the proof an L 1 —L°° decay estimate for solutions of the linear diffusion- 
advection 

dT 

— + V .\7T = AT (19) 

that is independent of the advection strength. Equation (|19|) is considered in the infinite strip 
D = R x [ — 7rZ, 7rZ] with the 27r/-periodic boundary conditions in y direction. 

Lemma 3.1 There exists a constant C > so that the solution of 

^ + v .VtP = AtP (20) 
at 

V>(0,x) = Vo(x) > 0, x£R 2 

with the 2nl -periodic boundary condition in y and a flow v that is 2irl -periodic, sufficiently regular 
and divergence-free: V • v = 0, satisfies 



L°° 



(D)<Cn 2 (t)U \\ L ^D), (21) 



10 



where D = M x x [— irl,Trl]y. Here n(t) is the unique solution of 

4n 4 (t) C x 



(22) 



l + 4n 3 (t)l 3 Pt, 
and the constants C,C\ do not depend on v. 

Remark. Note that ()22|) implies that for t > I , we have n{t) ~ . Hence, solution decays 
as the solution of the one-dimensional problem after the time it takes the diffusion to feel the 
boundary. 

Proof. We multiply (120H by ip and integrate over the domain D to obtain 

lj t u\\i = -\m\i (23) 

Here and below in the proof of this Lemma, || • || p denotes the norm in L P (D). 

We now prove the following version of the Nash inequality j^Tj for a strip of width I: 

'™' 5 £ c ^fMm <24) 



The proof of (|24j) is similar to that of the usual Nash inequality. We represent ijj in terms of its 
Fourier series-integral: 

where 

Mk) = \ [ e- ikx - in y/ l ^(x,y)dxdy. 
1 Jd 

Therefore we have \ifi n (k)\ < — \\iPWl~ 1 -- The Plancherel formula becomes 



D 



z,y)\ 2 dxdy = [ e in yl l - im y/ l+ikx - i P x ^ n {k)^ m { V ) dMpdxdy 

= iJ2[ l&(*)l 5 



(2tt)' 

dk 



and similarly 



jf \V^(x,y)\ 2 dxdy = lJ2 J m (* 2 + |^ n (fc)| 2 



(27T) 



2 ' 



Let p > be a positive number to be chosen later. Then using the Plancherel formula we may 
write 



2 2 = I + II, 



where 



\n\<pl J \ k \- p 
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The rest may be bounded by 



ne 

Therefore we have for all p > 0: 

i^< Cp(I ; +1> itf+ii^ 

We choose p so that 



< 



and obtain 

^ll^llfllV^Hf + C/ 2 /3||V^||f.-, 2 /3 



This is a quadratic inequality ax 2 +bx — c > with x = ||\7-i/>|| 2 / , a = l 2 ^ 3 \\ip\\i > & 
and c = HV'lli/C an d hence 



-6 + V6 2 + 4ac 2c c 

2a 6 + Vb 2 + 4oc _ Vb 2 + Aac 

This implies that 



l|V^||f >C M H^-+4^ M f||^|| 2 



8/3 \ - 1 / 2 



and therefore 



8/3 _,„ \ _3/2 /II./.II4 \-l 



> 



ci 2 H\\l 



Hence (|24|) indeed holds. 

We insert ()24|) into the inequality (|23j) and using the conservation of the L 1 -norm of i/j (recall 
that the initial data is non-negative) obtain 

" <- - M ° P 'g ,„ ■ (-) 



* " Uo\\f + l 3 \ 



y 0||l||^|l2 
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Integrating (|25jl in time we have 

lb n \\i Z 3 IU/)nlU 1 



where z(t) = \\ip(t)\\2/\\tpo\\i, an d thus 



4z 3 (t) 



4z 4 (i) 1 , s 

< 77727- ( 2 6) 



i + 4/ 3 z 3 (t) - cm 

The function on the left side of (|26|) is monotonically increasing and hence we have 

|2 < n(t)|hfo||i, (27) 



where n(t) is the solution of 1)22(1 . 

Let us denote by Vt the solution operator for (|20|): ^(i) = Vtipo- Then lj2T|) implies that 
II^IIli— >L2 < n(i). The adjoint operator T 7 ^ is the solution operator for 

^ - « • = (28) 
^(0,x) =^o(s), ieR d 

Note that the preceding estimates rely only on the skew adjointness of the convection operator 
v ■ V. Therefore we have the bound \\Vt W^—^l 2 — n {t) an d hence ||7 ? t||ji2-*£,oo < n(t) so that 

< n(t/2)U(t/2)\\ L 2 < n 2 (t/2)U \\ L i (29) 

and the proof of Lemma 13. II is complete. □ 

A very similar argument leads to an estimate for the solution of (|19|) with the 27r/-periodic 
or zero Dirichlet boundary conditions in both x and y. We state this variant which we will need. 

Lemma 3.2 Consider equation \2U\) with the 2ttI -periodic boundary conditions in x and y. As- 
sume that the initial data ipo(x) is mean zero: J i/jo(x)dx = 0. Then there exists a constant 
C > such that 

U(t)\\ L ^(D) < Cn 2 (t)\\4, \\ L i {D) , 
where D = [0, 2ttI] x x [0, 27r/]y. Here n(t) is the unique solution of 

4n 4 (t) Ci_ 
l + 4n 3 (t)P ~ Pt, 

and the constants C, C\ do not depend on v. The same result holds for the zero Dirichlet boundary 
conditions, without the assumption that the initial data is mean zero. 
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3.2 Variation of temperature on streamlines and gradient bounds 

The next step of the proof is to estimate the time-space averages of the oscillations of the solution 
along the streamlines of u and show that they are small if the advection amplitude is sufficiently 
strong. This does not have to be true in general, but a bound of this type holds if the initial 
data is uniform on the streamlines - then one has only to show that no large oscillations appear 
at later times. First, we reduce the problem to such initial data with the help of the Nash 
inequality in Lemma 13. II The maximum principle implies that it suffices to prove Theorem 11.21 
for the initial data of the form Tq = 1 for \x\ < Lq and To = elsewhere. We split the initial 
data as 

T (x) = T (x)r ? (/ l (x)/5o) + T„(x)(l - ^WA))) = 6>i + 002 ■ (30) 
The small parameter 5q will be specified later. The cutoff function 17(h) satisfies 

< T](h) < 1 for all heR, rj(h) = for \h\ < 1, r)(h) = 1 for \h\ > 2. 

We split the solution 4>(t,x) of (fTH|) as a sum <j) = 4>i + 4>2- The functions <f>i^ satisfy ((TH|) with 
the initial data </>oi and </>02, respectively. The function ^oi is smooth and is constant (1 or 0) 
on the streamlines of the flow. In particular it is equal to zero in the whole "water-pipe" system 
of boundary layers around the separatrices. The function ^02 satisfies a bound 

\\<\>m\\v(p) <|{xefl = lx [-Til, irl] : \x\ < L , \h(x)\ < 25 }\ < C5 L ln(l/5 ). (31) 

Therefore, Lemma 13 . 1 1 implies that if Lq is sufficiently large and we require that 

L (ln(L /0) 

with an appropriate constant C, then the function 4>2 satisfies a uniform upper bound 

\\Mt = l 2 )\\L~ { D)<^. (33) 

Hence we choose 5q as in (|32j) and concentrate on the function <pi that is initially uniform along 
the streamlines of u. We drop the subscript one to simplify the notation wherever this causes 
no confusion. 

Next, we obtain some uniform estimates for solutions of ()18|) that are initially constant on 
the streamlines. 

Lemma 3.3 For any time t > we have 

[ [ \V<p\ 2 d^ds < I |0 o (x)| 2 dx. (34) 
Jo Jd Jd 

Assume in addition that the initial data (fio(x) for the equation hltty are constant on streamlines. 
Then 

I [ \u-Vcp\ 2 d*ds<CA- 2 t [ {A^dx + CA- 1 [ \<po\ 2 dx. (35) 
Jo Jd Jd Jd 
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Proof. Multiplying (|18jl by <p(t, x) and integrating we trivially obtain 



|<X*,x)| 2 dx + 



D 



\V(f)(s,x)\ 2 dxds 



JD 



b (x)\ 2 dx, 



D 



which implies ()34j) . Next, we multiply (|18() by u • and integrate over D to get 
/ \u-V(f>\ 2 dx = A' 1 [ u ■ V0(A0 - </) t )dx 

JD JD 

< — ^ / 4>fdx + ^ I \u-V(f>\ 2 dx- A^ 1 I V4> T VuV(f>dx. 
2A J D 2 J D J D 



Thus we get 



/ \u-V<j)\ 2 dx<CA- 2 [ tfdx + CA- 1 [ \V0\ 2 dx. 

JD JD JD 

From (|36|) and (|34|) we obtain that 

I I \u-V^{x,s)\ 2 dxds <CA~ 2 [ [ ^dxds + CA' 1 

Jo JD JO Jd 



\ 2 dx. 



(36) 



(37) 



Recall that initially 4>o(x,y) is constant on streamlines of u, so that u ■ V0o = 0. Since the 
function 4> t satisfies the same equation (|TH|) as cj), this implies that 

/ |^(t,x)| 2 dx+ / / |V<Ms,x)| 2 dx(fe = [ |^t(0,x)| 2 dx= / |A0 o | 2 dx. 
Jd Jo Jd Jd Jd 

Therefore, bounding from above the time derivative term in ()37(l . we arrive at (|35|) . □ 
Note that the initial data for the function cpi in (|3fl|) obeys an upper bound 



D 



\A<j) 01 \ 2 dx 

-L JO 



L 



-L JO 



A 



Mx) 

So 



dydx 



A/i(x) ,(h{x)\ |V/i(x)| 2 „(h{x 



+ 



*o 2 



dxdy < „ In 
"o 



Hence, according to (|33|). the total oscillation on streamlines is bounded by 

/ [ |n-V0(s,x)| 2 dxds <CA~ 2 t [ lAcbo^dx + CA- 1 [ |^ i| 2 dx 
Jo Jd Jd Jd 



(38) 



< 



CA~ 2 TLr 



fi 3 



In ( 1 ] +r'J- J L ll 



for t < t. 
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3.3 Boundary layer and cell-to-cell heat conduction 

Lemma 13.31 suggests that for most times, there should be very little temperature variation along 
the streamlines. Our next goal is make this statement more precise. Let us fix a time r to be 
chosen later. We will denote by h and 9 the coordinates inside each cell, with h being our usual 
stream function and 9 the orthogonal coordinate normalized by the condition | V0| = |V/i| along 
the cell boundary and increasing in the direction of the flow. Let 5 > be an arbitrary, small 
number to be chosen later. 

Lemma 3.4 There exists ho satisfying 25 > ho > 5 such that 
[ E su Pex,e 2 10(^o, Out)- <j>{h , 9 2 ,t)\ 2 dt < CS-^LqI In (?j (A- 1 t5 3 In (J-^j +1) (39) 

where 5q satisfies h32fy. As a consequence, given a small number 7 > 0, for all times except for 
a set 5 7 of Lebesgue measure at most jt we have 



£ su P fc jft X/io, 0i, t) - Hho,02, t)\ 2 < CS^A'^Lol In (J\ (A~% 3 In (±-\ + r" 1 ) 

cells V / 

Proof. It is easy to check that 

^ = \Vh\- 1 \V9\- 1 u-V^. 
Let us denote by S a the streamline h = a £ [5, 25] in a given cell. Then we have 

S up0 1)0 J0(Mi)-<MM2)| 2 < (7 lu-V^ —^— ) 2 

„ f 1 de f d9 I f . ^ J|2 d9 



(40) 



The last step follows from the estimate 



is, |V0||Vfc| - \S, 
in the tube 5 < h < 25. Integrating in h and in time and summing over all cells, we obtain 

f [ ^sup gid2 \<j)(h,9 1 ,t)-<p(h,9 2 ,t)\ 2 dhdt<Clln(j S ) [ [ \u ■ Vcp\ 2 dxdt. (41) 

JS cells ' ^ 6JJ ° JD 

Now ()39|) and (j40|) follow from (|38[) by an application of the mean value theorem. □ 

We see that for all times but a small set of "exceptional" times, the value of the temperature 
on the h = ho streamline is close to some constant in any cell. Our next goal is to establish 
a control on how different these constants can be for two neighboring cells. Consider two such 
cells, C- and C+, and let B be their common boundary. Let us choose the coordinates on these 
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two cells so that h = on B, h > on the right cell C+, h < on the left cell C_, and the 
angular coordinates 6L and 9 + are equal to zero in the mid-point of B. Denote by D± the region 
bounded by curves #-,+ = ±y an d streamlines /i = ±/tQ, where y ~ / is chosen so that D-t- is 
sufficiently far away from the corners of the cells (so that |V/i|, |V#±| > C on -D±). Also denote 
by the pieces of streamlines bounding D±. Let us define the temperature drop between C_ 
and C + as follows: 

\4> + — </>_| = max |0, mins + </> — max5_<^>, mms_4> — max<j + </>} ■ 

Note that if the time is not exceptional, Lemma 13.41 implies that maximum and minimum of <ft 
along h = ho streamline in any cell differ by at most Cd^A'H / 5)^- 1 L {A- 1 5^ \n(l/S ) + 
r -1 ). Now we are ready to state 



Lemma 3.5 For any r > 0, we have for D = R 

10+ - <P-\ 2 dt < csr 1 

cells 



-ttL IT I] 



\Vc/)\ 2 dxdt < CL 5. 



(42) 



JD 



Therefore, given 7 > 0, for all times with an exception of a set of Lebesgue measure at most jr, 
we have 



cells 



(43) 



Proof. The proof is straightforward and has already appeared in |2()j . We sketch it here for the 
sake of completeness. Clearly, 



ho 



-ho 



dh 



(h,9)dh 



> 



for any 9± = 9. Since in the region D± we have |V/i|, |V0±| ~ 1, integrating in the curvilinear 
coordinates we obtain 



|V^| 2 dx > 



D 



y rho 
-y J— ho 



d(f> 



dh 



d9dh >CS~ l l\ 



finishing the proof. □ 

Lemmas 13.41 and 13.51 allow to control how much the temperature changes from cell to cell. 
Our next lemma summarizes this in a way that will prove useful. 

Lemma 3.6 Assume that at a certain time t, the estimates and hold. Suppose that 
there exists a cell C such that 4>(ho,9) > [3 > at some point (ho, 9) € C. Then for at least N 
cells, we have 4>(ho,9) > (5/2 for any 9 inside these cells, where the number N can be estimated 
from below as follows: 



N > Cp 2 jr (5L + S-^-A-^LqSqH ]n(l/S)r ln(Z/<S ) + S^A^LqI ]n(l/S)) 1 



(44) 
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Proof Consider a cell C\ which is the closest to C and such that there exists a point (ho, 9) G C\ 
with tf>(ho, 9) < (3/2. Then we must have 

Y / \^+-^-\ + J2 su ^ 1 ,eM h o,^)-Hho,e 2 )\>/3/2> (45) 

cells cells 

where the sum is over all cells between C and C\. On the other hand, Lemmas 13 . 41 and 13 . 51 imply 
that 

\<t>+ - ^-l 2 + E su ^ u e 2 \</>(ho,Oi) - 4>(h , 9 2 )\ 2 (46) 

cells cells 

< C-f- 1 T- 1 (5L + 6- 1 A- 1 L lln(l/5)(A- 1 6 3 Tlii(l/6o) + 1)). 

Since N^ =1 a 2 n > fem=i a «) 2 . a combination of gEJ) and © yields (gU). □ 

Now we can explain the strategy of the proof of Theorem ll.2l in more detail. Assume that at 
a certain " good" (not exceptional in the sense of Lemmas 13.41 13. 5 I) time we have a sufficiently 
high temperature (ft = (3 in some cell on a streamline h = h§. Then an appropriate choice of 5 
(hence ho) and an application of Lemma 13.61 ensure that (ft > (3/2 in many cells. Assume that 
for a sufficiently large portion of times < r, the temperature at h = ho streamlines is high in 
many cells. Clearly, the interiors of the cells will heat up too, giving (for a suitable choice of 5 
ensuring N is large enough) a contradiction with the preservation of the L 1 norm of (ft. On the 
other hand, if the temperature on the streamline h = ho is low most of the time, we expect the 
solution inside the cells to be small and quenching to happen if the cells are sufficiently small - 
the last condition ensures that the interaction between the boundary and the interior of the cell 
happens on a short time scale. To make the above plan work, we need a good understanding of 
an auxiliary "cell heating" problem. The next section is devoted to this goal. 

3.4 An auxiliary cell heating problem 

Consider a domain £1 which is a sub-domain of a cell, defined by the condition \ h\ > \ho\, that is, 
Q is the interior of the closed streamline {h = ho}- We consider the following initial-boundary 
value problem on Q : 

w t + Au ■ Vw - Aw = (47) 
w(t,x) = a(t), x G dil 
w(x,0) = g(x). 

Here a(t) is some (smooth) function of time, independent of x. Let us derive a compact formula 
for the solution of (|47|) . Set v = w — o~(t). Then v satisfies the zero Dirichlet boundary conditions, 
and we have 

v t + Au ■ Vv - Av = -a'(t), u(x, 0) = g(x) - a(0). 

Let us denote by TC the differential operator — A + Au ■ V with the zero Dirichlet boundary 
conditions on <9f2 and by e _tw the positive semigroup generated by 7i. Using the Duhamel 
formula, we get 

v(t,x) =e~ m (g-a(0)l) - [ e-^^la' (s) ds. (48) 

Jo 
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Here the semigroup is applied to a function equal identically to one in all of f2, denoted by 1. 
Integrating by parts in (|48|). and recalling that w = v + cr(t), we obtain 



w(t, x) = e 



m 5 (x)+ f a(s)He-^ n lds. (49) 



It will be convenient to use probabilistic interpretation of (|49|) . Recall that we can represent 
the semigroup action in the following way (see, e.g. |12j^: 

e-* w /(x) = E x [f(Xf (u,))] , 

where X* is the diffusion process corresponding to H starting at x, and the expectation E x 
is taken over all paths that start at x and never leave Q up to time t. The latter restriction 
corresponds to the zero Dirichlet boundary condition on d£l. Note that in the case / = 1, the 
function Q(t,x) = e~ tH \ coincides with probability that the stochastic process Xf{oj) never 
leaves before time t : 

p -m 1 = E x [!(x t x (o;))] = P(Xf{oj) G £2, Vs<(). 
In terms of the function Q(t,x) expression (|49|) becomes 

w(t,x) = e" m ^(x) - / cr(s)5tQ(t-s,x)ds (50) 

JO 

In order to control the solution of an auxiliary cell problem, we need to control the properties 
of Q(t,x), uniformly in the flow amplitude A - this is akin to Lemma \'A.1[ The first result we 
need in this direction is the following lemma, showing that the exit probability is bounded from 
above by a constant independent of A. 

Lemma 3.7 For any t satisfying < t < C\l 2 we have 

/ Q(t,x)dx > C 2 l 2 > 0, (51) 
Jn 

where C 2 is a universal constant depending only on C\ but not on A or I. 
Proof. Let us rewrite the generator 7i in the natural coordinates h, 9 : 

« = W£ + \Vh\^ + Ahl + (AO - 4V«||V»|)|. 

Note that Ah = —2l~ 2 h - this property is by no means crucial for our analysis but it does 
simplify some computations. The diffusion process Xf corresponding to Tt, written in the 
(h, #)-coordinates, is given by 

dX? = V2\Vh\dB^ ] - 2r 2 X?dt (52) 
dXf = V2\V9\dB<[ 2) + (AO - A\Vh\\V9\)dt, 
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where the values of all functions are taken at a point (Xj},Xf), while and P t (2) are inde- 
pendent one-dimensional Brownian motions. Clearly, Q(t,x) > P(X« G [0, Z], V s < i). It is 
not difficult to see that the exit probabilities of X^ are majorized by the exit probabilities of the 
Ornstein-Uhlenbeck process where the factor \Vh\ in ((o^l) is dropped. Indeed, let us introduce 

a(t,u) = V2\Vh(X?(u;),X e t (u J ))\. 

Multiplying (|52|) with e 21 2 * and integrating leads to 



e 2l-H x K*) = h(x) + f e ^- 2 s a{s ^ )dB {l) 

Jo 



Making a random time change (see, e.g., |29j). we find that the integral above has the same 
distribution as the Brownian motion B®, t ^ , where 

p(t,u) = / e 41 2s \a(s,uj)\ 2 ds. 
Jo 

Since |a(s,w)| 2 < 2, we have 

p(t, W )<p(t)^i/ 2 (e 4/ - 2t -l), 

and therefore 

Q(i,x) > P(B h p $e- 2r2s e [0,1], V s < t). (53) 

Let us remark that the expression on the right hand side is exactly the exit probability for 
the Ornstein-Uhlenbeck process mentioned above. The claim of the lemma now follows from a 
simple rescaling r = t/l 2 , Y T = l~ 1 Xp T . □ 

Next, we need some information on the behavior of dtQ(t,x). 

Lemma 3.8 For any t satisfying C\l 2 > t > 0, we have 

[ d t Q(t,x)dx < -C 3 , 
Jn 

with a constant C3 > that depends on C\ but is independent of A or I. Moreover, J n dtQ(t, x)dx 
is monotonically increasing in time. 

Proof. According to the previous lemma, 

/ Q(t,x)dx = \\Q(t,x)\\ L i > C 2 l 2 > 
Jn 

for any t < C\l 2 . The Cauchy-Schwartz inequality then implies a lower bound for the 1? norm: 
||Q(t, x) 11^2 > C2L Recall that the function Q(t,x) solves 

^ = HQ(t,x) 
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with the zero Dirichlet boundary conditions and the initial condition Q(0,x) = 1. Then 

^||Q||i 2 = -||VQ||i 2 <-r 2 ||Q||i 2 (54) 

by the Poincare inequality. Thus we have J QQtdx < — C| where C2 is independent of A and I. 
It is clear that < Q(t, x) < 1, and also that dtQ(t,x) < since Q(t,x) is the probability that 
the diffusion process starting at x does not exit £1 before time t. This implies the first statement 
of the lemma. 

To prove the second statement, we use again the monotonic decay of Q(i,x) in time. Fixing 
some time t > and integrating by parts (using the fact that the flow u is tangent to <9f2) we 
find that 

/ d t Q(t,x.)dx = -^-ds. 
Jn Jan on 

Since the function Q vanishes on the boundary and decays (in time) inside, we see that J n dtQ dx 
is monotonically increasing in time. □ 

Corollary 3.9 The solution w(t,x.) of J^7| ) with any nonnegative boundary and initial data 
satisfies 

/ w(t,x)dx>C / a(s)ds, (55) 
Jn Jo 

where C is a positive constant depending on t but not on A. If t < C\l 2 , the constant C can be 
chosen independent of I. 

Proof. The corollary follows immediately from the previous lemma and the representation 1)50(1 . 
□ 

Since we have no control over what happens at a small set of exceptional times, we need an 
estimate from above on how much the L 1 norm of the solution can change if the boundary data 
is close to one for a short time. The second part of Lemma 13 . 81 implies that it is sufficient to look 
at the situation when this hot period occurs at the end of the interval [0,t]. More precisely, if 
we replace a(t) by its monotonically increasing rearrangement then the L 1 -norm of the solution 
increases. The following lemma will be useful in such scenario. 

Lemma 3.10 For a time t satisfying < t < I 2 , we have 

j (1 - Q(t, x)) dx < Cl 2 {tr 2 fl 2 In (j) . (56) 

Proof. Note that 

/(1-Q(i,x))dx= f P(3s<t: Xf(co) £Q)dx. 
Jn Jn 

It follows from the proof of Lemma 13.71 that the probability on the right hand side is majorized 

by 

P(3 s < t : B^e- 2l ~ 2s $ [0, 1]) < P(3 s < t : $ [0, 1]), (57) 
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where p(t) = \l 2 (e il t — 1). Take a constant C so that p(t) < Ct for any t as in the statement 
of the lemma. Let t = tl~ 2 be the rescaled time. Then the probability on the right hand side of 
(|57|) does not exceed the probability that B® leaves the interval C^ 1 / 2 [— h/l, 1 — h/l] before the 
time t. Using the reflection principle, we can estimate this probability from above by one, if h/l 
or 1 — h/l are less than Vt, and by C n (tl/h) n otherwise (where n > is arbitrary). Integrating 
and taking into account the Jacobian of transformation between h, 9 and x, we arrive at 




which is nothing but (|56|) . □ 



3.5 Proof of Theorem fOl 

Now we are ready to complete the proof of Theorem 11.21 In this proof, we will work with the 
time scale r = CqI 2 for a sufficiently large universal constant Co- It means that in all estimates 
of the previous sections that we are going to use, we take r = Col 2 . This ensures, in particular, 
applicability of Lemma 13.81 and Corollary 13.91 If we show that ^(i,x) < 9' /2 = #o/(2e) for some 
t < t, this will be sufficient for quenching since we assume that r < r c = M _1 (small cell size). 
Proof. We have a set 5 7 of exceptional times of size at most jt such that on the complement 
of S 7 the estimates (|40|). (j4"3]l are valid. The value of the constant 7 will be chosen later and 
will be independent of A and Assume first that for all t G except for a set Sb (of "bad" 
times) of size at most 7T we have cp(ho,9) < 7 for all cells and all 9. We are going to show 
that if 7 is chosen sufficiently small (with the choice being uniform in A, I), then we must have 
quenching. Let fi n be the subset of a cell C n enclosed by \h\ = \ho\. Then for any x e Q n we 
have </>(£, x) < <fi(t,x) where cp(t,x) satisfies 

4> t - A(f> + Au ■ = 

with the initial data </>(x, 0) = 1 in O n , and the boundary data given by 

Ms x)l = 1 1 ' S £ 57 U Sb ' 1 
' I7, otherwise 

This bound on <p follows from the representation formula (|50j) and the above assumption on the 
behavior of 4>(ho,6). By linearity, inside each region Q n , we have 4>(t,x) = <f>i(t, x) + 02(t,x), 
where 4>i satisfies the zero Dirichlet boundary conditions, while 4>2 has zero initial data. Let 
t\ = C\l 2 , where C\ is a universal constant which will be chosen below. A simple argument 
using (|54|) shows that 

r 2 \\Mti,-)\\mn n) <r 1 ||$ 1 (t 1 ,oili2 (nri) (58) 

uniformly in A, / as C\ — > oo. 

The function 02 can be estimated in the L 1 -norm using (cf. (|50^ 

ft 



4> 2 (t,x) = — a(s)d t Q(t - s,x)ds 
Jo 
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and the second statement in Lemma 13.81 as 

\\<h(ti, •) || v- (Q„) <- / cr(s) / d tl Q(h- s,x)dxds 



o 



/•tl /• /"El — p 7 U0(,| /■ 

<— / / d tl Q(ti — s,x)dxds — 7 / / d t Q{t\ — s,x)dxds 

< / [l-Q(|5 , 7 U5 b |,x)]dx + 7 f [1-Q(ti,x)]dx 

(59) 



< C/ 2 



K:)) 1/2 H^ ) -K'n 



by Lemma 13. 1UI 

It follows now from the parabolic maximum principle that for all times i > ii, the solution 
x) of the original linear advection-diffusion problem l|18|l in any cell C n satisfies 

0(t,x) < 0*(t,x) +0 2 (*,x). 

The function 02 has been estimated in (|31|) and (|33|) , The function x) solves (|47j) for t>t\ 
with the 27r/-periodic boundary conditions in x, y and at time t = t\ = C\l 2 is given by 



**(ti,x) = { f ljX) ' X ^ 

^ 1 X t Cfj 



x G C n \ ^ra> 

Lemma 13.21 allows us to choose a positive number 7 so that if 

ii<m,x)ii L1(c „)<7^ 

then 

+ i 2 ,*)\\ L - {Cn) < Cn 2 (/ 2 )||<f (^x)!!^ < e'jA. 

This is possible, as n 2 (/ 2 ) ~ /~ 2 (see Remark after Lemma I3.1JI . Note that 

||0*(ti,x)|| zl(Cn) < ll^iCii.xJiUx^ + n^Cti.xJiUi^ + c?*/^/*). 

Choosing Cx large enough we can use (|58|) to estimate the first term on the right hand side 
and make sure it does not exceed 7^ 2 /3. Next we choose 7 small enough so that (|59|) gives 
similar control of the second term. Finally, choose 5 small enough so that the last term is also 
sufficiently small. With this choice of C%, 7 and 5 (uniform in A, I) we have 

Ufa + Z 2 ,x)|| L ~ < U*(t! +/ 2 ,x)|| i o 0(Cn) +9' /W < 9' /2 

and thus quenching. 

It remains to consider the case when there exists a set of bad times Sb E [0, ti] of size at 
least 7T (with 7 universal constant determined above) such that for any t S Sb, there exists a 
cell C n such that for some point (ho, 6) £ C n , we have 4>(ho,9) > 7. We are going to show that 
this cannot be true if A is large enough, thus forcing the scenario which is considered above and 
leads to quenching. Indeed, Lemma 13.61 implies that there are at least 

N = C7 3 r (5L + 5~ l A~ 2 tI \n{l/5) In (l/6 ) + A~ l L Q l\n(l / 5)) _1 (60) 
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cells such that (ft(ho,6) > 7 /2 for any 9 in these cells. For each cell, let a(t) = ming(ft(ho,9,t). 
On each cell, solve the initial-boundary value problem (|47|) with g = 0, and denote the solution 
4>(t,x). By the parabolic maximum principle, we have that (ft(t,x) > cft(t,x) for \h\ > |/io|- 
Applying Corollary 13,91 we obtain 

/ (j)(r, x)dx>CY^ / (») ds = C a(s) (is > C 7 2 riV, (61) 

^ ceiis ^ " ce/Zs 

where iV is given by (|6U|) . We claim that if <5 = <5(Lo) is chosen to be sufficiently small and 
A = A(Lq) sufficiently large then (ffiTj) leads to 

J(r, x) dx » ZL , (62) 



obtaining a contradiction since the L 1 norm of (ft is preserved by evolution. Indeed, since 7 and 
t/Z 2 are fixed constants, it suffices to choose S so that 



T 2 



(63) 



and then choose A so that 



Recalling the formula (|32"|) for 5q, we discover that A > C(1)Lq\ti Lq is sufficient to satisfy (|64j) . 
H63|). completing the proof of Theorem ll.2l □ 
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